Seismic inference of differential rotation in Procyon A 
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Abstract. The differential rotation of the F5V-IV star Procyon A is computed for a class of models which are consistent with 
recent astrometric and asteroseismic data. The rotation pattern is determined by solving the Reynolds equation for motion, 
including the convective energy transport, where the latter is anisotropic owing to the Coriolis force action which produces a 
horizontal temperature gradient between the poles and the equator. All the models show a decrease of the rotation rate with 
increasing radius and solar-type isorotation surfaces with the equator rotating faster than the poles, the horizontal rotational 
shear being much smaller for models with a less extended convective envelope. The meridional flow circulation can be either 
clockwise or counter-clockwise, and in some cases a double latitudinal cell appears. The rotational splittings are calculated for 
low degree p-modes with I = 1, m = 1 and I = 2, m = 1, 2, and it is shown that, for modes with m = 1, the stronger is 
the horizontal differential rotation shear the weaker the effect on the average rotational splitting expected, whilst the opposite 
happens for the mode with m = 2. On the basis of the present study, a resolution of 10 nHz in individual oscillation frequencies 
seems to be necessary to test the different dynamical behaviour of the proposed models, that appears barely achievable even in 
the forthcoming space missions. However, the average over several observed splittings could produce the required accuracy. 
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1. Introduction 

One of the brightest stars in the sky (my = 0.34) is Procyon 
A (a CMi), hereinafter simply called Procyon, an extensively 
studied object whose evolutionary phase is believed to be that 
of a subgiant star with a convective core, a radiative enve- 
lope and a thin surface convective zone, where p-modes are 
presumably excited (Christensen-Dalsgaard & Frandsen 1983; 
Houdeketal. 1999). 

For these reasons Procyon is a favorite target for astero- 
seismic observations. The first clear evidence for the presence 
of oscillations was obtained by Martic et al. (1999), but when 
it was observed in 2004 by the MOST space mission, the re- 
ported precision from the photometric measurements (15 ppm) 
was not sufficient for the detection of p-modes (Mattews et al. 
2004). On the other hand, three different teams have then re- 
ported the identification of p-modes on Procyon (Martic et al. 
2004; Eggenberger et al. 2004; Claudi et al. 2005), leading to 
the determination of the large frequency separation in the range 
54 /iHz - 56 //Hz. It is thus conceivable that future observation 
campaigns from ground or space will further refine the deter- 
mination of the basic asteroseismic parameters and, in addition, 
will be able to measure a possible rotational splitting in this 
slow rotator. 
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The present work is addressed to the determination of the 
differential rotation for a class of models of Procyon, consistent 
with the presently available astrometric and asteroseismic data, 
and the related rotational splittings. As the evolutionary mod- 
els of Procyon show that the mass in the surface convection 
zone is very small, ranging from ~ 10~ 7 M@ to ~ 10~ 4 M Q 
during the evolution, it is not clear which type (if any) of differ- 
ential rotation is present. In particular, the question we would 
like to investigate is whether a significant amount of latitudi- 
nal shear is present, so that a solar type of differential rotation, 
~ 0(r, 8), is generated, or instead O ~ fi(r), as is often as- 
sumed in asteroseismic investigations (Di Mauro et al. 2000). 

Recent advances in mean-field theory of solar differential 
rotation have proven to be very successful in explaining inter- 
nal solar rotation. The solar rotation pattern can be explained 
in terms of angular momentum transport by the gas flow in the 
convection zone. Assuming that Reynolds stresses are the only 
relevant mechanism of angular momentum transport, Kiiker 
et al. (1993) found a remarkably good agreement between 
the results of their model and the "observed" rotation, as de- 
duced from the inversion of rotational splittings. Kitchatinov 
& Rudiger (1995) and Kiiker & Stix (2001) presented a refined 
model that included both the meridional flow and the convec- 
tive heat transport, and showed that the flows driven by differ- 
ential rotation and a small horizontal temperature gradient, ow- 
ing to the anisotropic convective heat transport, have opposite 



directions, and roughly balance each other in the solar convec- 
tion zone. This analysis was extended to a generic class of zero 
age main sequence F stars of mass M = 1.2 Mq in which the 
convection zone was about 20% of the stellar radius (Kiiker & 
Riidiger 2005a) leading to the conclusion that a solar-type rota- 
tion is the most common rotational law, at least for moderately 
rotating solar type stars. 

In the case of Procyon one instead deals with a M ~ 
1.42 Mq — 1.5 Mq evolved F star which presumably is end- 
ing or just ended the main sequence phase, with a convection 
zone much more shallow than that of the cases discussed in 
Kiiker & Riidiger (2005a). 

2. The evolutionary models 

The starting step of our investigation was to construct evolu- 
tionary models of Procyon by means of the GARSOM code 
(Schlattl et al. 1997, 1999; Bonanno et al. 2002) whose struc- 
ture provides the basic stratification for solving the Reynolds 
equation (Kiiker & Stix 2001). 

In the present calculations the OPAL opacities (Iglesias 
& Rogers 1996) implemented in the low-temperature regime 
by Alexander & Fergusson (1994) tables have been used. As 
far as the equation of state is concerned, we used either the 
OPAL-EOS 1996 (Rogers et al. 1996) or the MHD (Hummer 
& Mihalas 1988; Mihalas et al. 1988; Diippen et al. 1988) ta- 
bles. The OPAL-EOS 1996 tables have been updated by a rel- 
ativistic treatment of the electrons and by improving the ac- 
tivity expansion method for repulsive interactions (Rogers & 
Nayfonov 2002). The nuclear reaction rates have been taken 
either from Bahcall et al. (1995) or Adelberger et al. (1998). 
The convection has been treated by the classical mixing length 
theory. Microscopic diffusion of hydrogen, helium and all the 
major metals has been taken into account. The outer boundary 
conditions have been determined by assuming an Eddington 
gray atmosphere. 

The most recent determinations of the astrometric param- 
eters of the visual binary orbit of Procyon have been obtained 
by means of the Multichannel Astrometric Photometer obser- 
vations during the period 1986-2004. The inferred mass of the 
primary component is found to be M = 1.43 ± 0.034 Mq 
(Gatewood & Han 2006), whilst a previous determination by 
Girard at al. (2000) reported M = 1.497 ± 0.037 M®, with a 
parallax n = 283. 2 ±1.5 mas. Although the two measurements 
are consistent within the error bars, if only the Wide Field and 
Planetary Camera 2 observations are used in combination with 
the Hipparcos parallax, II = 285.93 ± 0.88 mas, the resulting 
mass is M = 1.42 ± 0.04 Mq, much closer to the Gatewood 
& Han (2006) value. Here we shall assume the conservative 
range 1.42 Mq - 1.50 Mq to constrain the mass of Procyon in 
our models, although the average value of these measurements 
would give a mass M = 1.45 ± 0.02 M . 

The luminosity has been estimated from the mean visual 
magnitude obtained by Allende Prieto et al. (2001) and bolo- 
metric correction by Flower (1996), so that log L/Lq = 
0.84 ± 0.02. Recent measurements of the angular diameter by 
Kervella et al. (2004) lead to R/Rq = 2.048±0.025, while the 
parallax of Girard et al. (2000) gives R/Rq = 2.067 ± 0.028. 



Ihe effective temperature has been calculated by ruhrmann et 
al. (1997) who found logT c = 3.815 ± 0.006. 



The metallicity of Procyon, as derived from spectroscopic 
measurements, results to be close to the solar one (Takeda et 
al. 1996; Kato et al. 1996), and for simplicity we adopted the 
same (standard) metallicity mixture of the Sun (Grevesse et 
al. 1996) so that at the surface (Z/X) ph = 0.0245. The in- 
clusion of metal diffusion in evolutionary models of Procyon 
leads to an almost complete depletion of the metal content 
at the surface (Chaboyer et al. 1994). The reason is that the 
gravitational settling is expected to take place below the con- 
vection zone which is only ~ 10~ 4 Mq thick. In order to 
avoid this problem, we have considered a mass loss accord- 
ing with the semi-empirical Reimers (1975) prescription, M = 
—4 x 10 -13 rj (L/gR), where the surface gravity, g, luminos- 
ity, L, and radius, R, are expressed in solar units, and r/ is an 
efficiency parameter, whose value depends on the evolution- 
ary phase of the star from the main sequence to the red gi- 
ant branch; here a value close to 0.1 results to be adequate. 
The rotational velocity, as estimated by Allende Prieto et al. 
(2002), is v sin i = 2.7 ± l.Okms -1 . If we assume that the ro- 
tational axis of Procyon is perpendicular to the plane of visual 
orbit (i = 31.1° ± 0.6°, Girard et al. 2000), the surface rota- 
tional velocity should be 5.2 ± 1.9 km s -1 , which corresponds 
to a rotational period of about 20 d. On other the hand, in the 
extreme case of i = 90° the rotational period is about 33 d. 
Present asteroseismic data are not accurate enough to put strong 
constraints on the computed models, although several groups 
have succeeded in measuring the large frequency separation 
Av . Martic et al. (2004) found Av Q = 53.6 ± 0.5 //Hz, while 
Eggenberger et al. (2004) A^o = 55.5 ± 0.5 fiHz. More recent 
measurements, obtained with the high resolution spectrograph 
SARG operating with the 3.5 m Italian Galileo Telescope at 
Canary Islands, report instead A^o = 56.2±0.5//Hz(Claudiet 
al. 2005). We thus constructed evolutionary models of Procyon 
constrained by 54 /iHz < A^o < 56 /iHz, and computed 
four representative models with masses in the range 1.42 Mq 
-1.5 Mq. In all the models, the the mixing length parameter, 
initial helium and metal abundances, diffusion of the elements, 
and the mass loss parameter ?; have been adjusted in order to 
reproduce the observed luminosity and effective temperature 
within a box delimiting their determination uncertainties, with 
the above mentioned constraint on Ai/q. As far as the variation 
of the mass loss parameter rj is concerned, an exploration in the 
range 0.05 < T] < 0.15, that is still consistent with the aster- 
oseismic requirements of matching the value of A^o, does not 
produce sensible effects on the stratification of the convection 
zone and by consequence on its dynamical behaviour. 



The basic characteristics of the four representative models 
are summarized in Tabled while the evolutionary tracks of the 
models are presented in Fig.^ Note that the Model 2 represents 
the model with a mass corresponding to the average mass of the 
most recent determination, as already discussed. 



lable 1. Characteristic quantities or the tour Frocyon models. Ihe suffixes ph and cz indicate top (photosphere) and bottom or 
convective envelope, respectively. Model 1 includes diffusion of H and He only, while Models 2, 3 and 4 include diffusion of all 
the elements. 
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Fig. 1. The evolutionary tracks of the models summarized in 
Tabled The solid line is for Model 1, the dashed for Model 2, 
the dot-dashed for Model 3, and the dot-dot-dashed for Model 
4. The box indicates the uncertainties in the observation deter- 
minations of luminosity and effective temperature. 

3. Differential rotation and heat transport in the 
outer convection zone 

The solar differential rotation was successfully reproduced by 
means of a mean-field model based on the Reynolds stress 
transport of angular momentum (Kiiker et al. 1993), improved 
by the inclusion of meridional flow and turbulent heat trans- 
port (Kitchatinov & Riidiger 1995). The inclusion of the above 
mentioned ingredients turned out to be crucial to obtain results 
in agreement with the isorotation surfaces deduced from helio- 
seismology, as the Taylor-Proudman theorem predicts cylindri- 
cal surfaces of isorotation in an isothermal fluid at high Taylor 
numbers. Since the Taylor number assumes values as large as 
10 s in the lower part of the solar convection zone, owing to its 
inviscid character, the state of rotation in cylinders would have 
been unavoidable (Riidiger et al. 1998). The model has been 
applied to several types of stars, either fully convective or with 
outer convection zones (Kitchatinov & Riidiger 1999; Kiiker & 
Stix 2001 ; Kiiker & Riidiger 2005a,b). 

A relevant quantity to understand the role of the turbu- 
lent transport of the angular momentum is the Coriolis num- 
ber Q* = 2rf2, where £1 is the stellar angular velocity and r 
the characteristic turnover time of convective eddies (Riidiger 
1989). The radial dependence of Q* in the outer layers of the 
star for a representative model of Procyon is shown in Fig.|2] 
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Fig. 2. Coriolis number as a function of the stellar fractional 
radius for a star with Al — 1.45 Mq in the subgiant phase for 
a rotation period of 20 d. 



Since the inhomogeneous L-effect depends on Q*^ 1 while the 
anisotropic L-effect depends on Q* _3 , it is then clear that they 
are both essential to correctly describe the turbulent transport 
of angular momentum (Kitchatinov & Riidiger 1995, 1999). 

3. 1. Reference model 

Since the influence of rotation and large-scale meridional flow 
on the structure of the star is very small we treat them as pertur- 
bations of an unperturbed, non-rotating, spherically symmetric 
reference state, that is a stellar convection zone in hydrostatic 
equilibrium, as computed by the stellar evolution code above 
described. In order to compute the differential rotation, we sim- 
plify the reference state by assuming a perfect gas with con- 
stant specific entropy and total heat flux throughout the con- 
vection zone. The differential equations for temperature and 
gravity (Poisson equation for a spherically symmetric gravita- 
tional potential with g(r) = — d<I>/dr) are those described in 
Kitchatinov & Riidiger (1995) for an analogous problem: 

dT = _s(r) 

dr c p 

^) = _ 2 ^) +47r Gp (2) 
dr r 



with the adiabatic relationship: 



T 



1/(7-1) 



(3) 



where c p — 7\L/V a d is the specific heat at constant pressure, 
with 1Z the gas constant and V a d the adiabatic gradient, g(r) 
the gravity, 7 = (1 — V a d) _1 the adiabatic exponent, and 
the suffix e indicates the values of the quantities at the outer 
boundary. The solution is determined by choosing a reference 
point from the model produced by the stellar evolution code, 
which provides the initial values necessary for the integration. 
Equations are then integrated outwards and inwards from the 
starting point at which the values from the reference model 
have been chosen, the inward integration being terminated at 
the point where the radiative flux reaches the total heat flux, 
which coincides with the bottom of the convection zone, and 
the outward integration being terminated at the stellar radius 
as determined by the reference model. Thus Eqs. Q and (|2jl 
constitute an initial value problem. To adjust the depth of the 
convection zone to the value of the reference model, we take 
the adiabatic gradient as input parameter. For fully ionized 
monatomic gas, V a d = 0.4. In the outer parts of the convection 
zones, however, it is smaller owing to the partial ionization. The 
stratifications found by the stellar evolution code for the mod- 
els of Procyon we considered show values of V a d as low as 
0. 1 in the outer parts of the convection zone. Since we assumed 
it constant with radius in our simple model, the values of V a d 
that give the best fit to the evolutionary models lie around 0.35. 
We thus produce a convection zone with adiabatic stratification 
that closely resembles the one reproduced by the evolutionary 
models, that is used as the unperturbed, reference state for the 
computation of the differential rotation. In spite of the simplifi- 
cations used, the deviations of our adiabatic stratification from 
that of the evolutionary model are within a few percent in the 
relevant quantities everywhere in the convection zone. 



3.2. Heat transport 

Following the procedure described in Riidiger & Hollerbach 
(2004), the mean-field equation of convective heat transport is 
given by the following expression: 



F rad ) + pc p u ■ (3, 



(4) 



where s is the mean specific entropy. The transport terms on 
the RHS of Eq. © are the convective heat flux, F l conv = 
pc p {v!{F') = pCpXijPj, the radiative heat flux, i^ rad = 
— lQaT 3 \/iT/3np, and the heat advection by meridional flow, 
with u the mean velocity field, where n is the opacity, a the 
Stefan-Boltzmann constant, and f3 = g/c p — VT = —V5T 
the super-adiabatic gradient. As the reference state has con- 
stant entropy sq, we define the perturbation as 5s = s — sq, 
and, since all derivatives of sq vanish, we can rewrite Eq. @ in 
terms of 5 s: 



F™ d )+pc p u-(3. 



(5) 



On applying the mass conservation law, V • (pu) = 0, re- 
placing the temperature perturbation with the entropy pertur- 
bation, 5s = Cp6T/T, we can rewrite the last term in Eq. (0 
as pc p u ■ (3 = — V • (pTuSs). Since we find in some cases that 
the heat transport by the meridional flow exceeds that of the 
convective flux, we introduce a horizontally-averaged entropy 
perturbation: 



5s\ = 



1 



5s sin 9d9 



thus splitting 5s into the two terms 5s~i(r) and 5s2(r,9). We 
then drop out 5s~i from the advection term to prevent convec- 
tive instabilities in our system, and obtain the equation of heat 
transport in its final form: 



rirad 



pTu5s 2 ) 



(6) 



3.3. Equation of motion 



The equation of motion for the mean velocity field u is (cf. 
Riidiger 1989): 



P 



8u 
~8t 



+ (ft ■ V)u 



-V • pQ - VP + pg 



(7) 



where Q is the one-point correlation tensor of fluctuating ve- 
locities Qjj = (u'ju'j). On assuming axial symmetry, the az- 
imuthal component of Eq. |7} can be written as: 



pr 2 sin 2 9— + V • t = 
at 



(8) 



where t = r sin 8(pr sin 6flu m +p(u' ( pu')). The two remaining 
components of Eq. describe the meridional circulation u m , 
and can be represented by a single scalar function, the stream 
function A: 



u = 



1 



dA 



1 dA 



pr 2 sin 9 89 ' pr sin 9 dr 







(9) 



On taking the azimuthal component of the curl of Eq. @ we 
obtain: 



8uj 
~dt 



1 



V x -VOQ) 



— (Vpx Vp), 



r sm 



dz 



(10) 



where ui = (V x u)<j, is the azimuthal component of the mean 
vorticity, d/dz = cos 9 ■ d/dr — sin6*/r • 8/ 89 represents the 
gradient along the axis of rotation, and dots indicate second- 
order terms of minor importance. As the reference state is in 
hydrostatic equilibrium, VP = — gp, we can rewrite the baro- 
clinic term on the RHS of Eq. (II Oi in terms of specific entropy: 



1 



-j(Vp x VP), 
P 



(Vs x VP). 

c p p 



^ (ID 

rc, 



8s 
89 



On neglecting the second-order terms, we have three terms on 
the RHS of Eq. (II 01 : the first, which contains the Reynolds 
stresses, is purely viscous, while the second and third drive the 
meridional flow. Close to the poles, the gradient of the rotation 



rate along the rotation axis dominates, as the horizontal entropy 
gradient vanishes there. At low latitudes, on the other hand, 
dft/dz vanishes and the baroclinic term dominates. As both 
terms can be either positive or negative, they can in principle 
enhance or cancel each other. For solar-type stars it has been 
found that both terms balance each other in a wide range of 
rotation rates, leading to a slow meridional flow and a strong 
differential rotation (Kiiker & Riidiger 2005). 

3.4. Transport coefficients 

Instead of the standard mixing-length theory we use the expres- 
sions derived by Kitchatinov & Riidiger (1999) for the convec- 
tive heat transport and Reynolds stresses. The heat conductiv- 
ity tensor, Xij< can be written as the product of a scalar dif- 
fusion coefficient, Xu an d a dimensionless tensor, $y, so that 
Xij = Xt$ij, where xt = c x T c ga 2 H 2 {/3 r )/4T, with r c the 
turbulent eddy correlation time, a the classical mixing-length 
parameter, H p the pressure scale height, c x a dimensionless co- 
efficient for heat transport, and the triangular brackets averages 
over 6. The stress tensor has the general form %j = —pQij, 
with the correlation tensor: 



= -M 



ijkl 



dn k 

dx. 



(12) 



The first term on the RHS vanishes for rigid rotation and repre- 
sents an anisotropic viscosity; the second one, the L-effect, is 
non-zero for rigid rotation, and therefore causes differential ro- 
tation unless balanced by the meridional flow. We write the vis- 
cosity tensor, Afijki — Vt^ijkU m a similar way as the convec- 
tive heat transport, namely v t = c v T c ga 2 H 2 (f3 r ) / AT , where 
c u is a dimensionless coefficient for viscous transport. The L- 
effect appears only in two components of the correlation tensor, 
and it can be written as follows: 



= Vt V sinOn, V = V° + V 1 sin 2 



COS 



v t HcosOVL : H = V 1 sm 2 9. 



(13) 
(14) 



The tensors <£>jj and ^ijki contain the anisotropy of the con- 
vective heat transport and the viscous part of the stress ten- 
sor, as computed by Kitchatinov et al. (1994). Together with 
the L-effect coefficients V° and V 1 (Kitchatinov & Riidiger 
1993, 2005) determine the transport of heat and angular mo- 
mentum by convection. Note that the updated L-effect coeffi- 
cients of Kitchatinov & Riidiger (2005) remove the discrep- 
ancy between theoretical and observed angular velocity dis- 
tribution in the outer 30, 000 km of the solar convection zone 
found in models which uses the older expressions. To adjust 
the model for the updated expressions for the L-effect, we now 
use c u — 0.27 with a — 2 rather than the values from Kiiker 
& Riidiger (2005a,b), who used c v = 0.15 with a = 5/3. With 
a Prandtl number, P r = Vt/Xt — 0-8 (Kitchainov et al. 1994), 
there is a very good agreement between the computed and ob- 
served rotation pattern for the Sun. The transport coefficients 
have been calibrated with the Sun, whose differential rotation 
is known from helioseismology. Numerical simulations on the 
effect of changes of transport coefficients show that the strength 
of horizontal rotational shear is quite insensitive to changes of 



parameters, while meridional flow is more sensitive leading to 
slight changes in the rotational pattern. Both the existence of 
the L-effect and the anisotropy of the convective heat trans- 
port have been confirmed by simulations (Kapyla et al. 2004; 
Riidiger et al. 2005). 

3.5. Boundary conditions and gravity darkening 

The numerical code we used to solve Eqs. 0, (|8), and (II Oi is 
the same as that used by Kiiker & Stix (2001). It uses a second- 
order finite difference scheme that conserves angular momen- 
tum and energy. The equations for the reference model have 
been solved by using a fourth-order Runge-Kutta scheme. In 
the equation of motion, we impose stress-free boundary condi- 
tions on both boundaries. In the heat transport equation, we fix 
the total luminosity on both boundaries, and assume that is in- 
dependent of latitude, so that L = Anr 2 F. This fixes the radial 
component of the entropy gradient, but not the entropy itself, 
so not precluding a horizontal entropy gradient. As we have 
written the heat transport equation as a conservation law, these 
boundary conditions ensure that a stationary state is eventually 
reached. 

Our code also includes the option of taking into account 
the gravity darkening by adding a rotation-dependent term to 
the heat flux equation: 



F = 



L 



47rr 2 



1 



2s 



3 2 a 1 
- cos 9 

2 2 



(15) 



where e = Q 2 r 3 /GM (Riidiger & Kiiker 2002). For a star 
like Procyon with a rotation period of 20 d, a value of e = 
1.8 x 10~ 4 is found, small value that cannot affect luminosity. 
In fact we did not find a significant difference between runs 
with and without gravity darkening. 

4. Rotational splittings 

The rotational splittings of oscillation frequencies depend on 
the radial and latitudinal angular velocity distribution in the in- 
terior of the star. It could be of interest, essentially in the light 
of possible future observations, to calculate the rotational split- 
tings for low-degree p-modes expected for our differential rota- 
tion models of Procyon. We use the following parameterization 
for the rotation law: 



n(r,9) = £ Qi(r) cos 21 9 



(16) 



so that the angular integration can be explicitly performed, and 
the rotational splittings are (Cuypers 1980): 



m * max f R 

with the kernels K n i m i (Christensen-Dalsgaard 1997): 

T^nlmi — 

Pr 2 {{{tr - a) 2 + (1(1 + 1) - 2* 2 - 3i- 2)&]Q 
+t(2i-l)Q /TO< _i}/J 



(17) 
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Fig. 3. The normalized angular velocity profiles for the four 
models with a rotation period P rot = 20 d. The four panels 
show respectively: top left: Model 1 ; top right: Model 2; bottom 
left: Model 3; bottom right: Model 4. The different line-styles 
denote different latitudes: solid: equator; dotted: 30°; dashed: 
60°; dash-dotted: poles. The normalization of the angular ve- 
locity is done with respect to an average angular velocity Cl as 
described in the text. 





Fig. 4. The meridional flow patterns corresponding to the ro- 
tation profiles in Fig. [5] The four panels show respectively: top 
left: Model 1; top right: Model 2; bottom left: Model 3; bottom 
right: Model 4. The solid lines of iso-velocity denote clockwise 
and the dashed ones counter-clockwise motions. For clarity, in 
the diagrams the depth of the convection zone has been en- 
larged by a factor of 2 for Models 1 and 2 and 4 for Models 3 
and 4. The meridional motion maximum speeds at the stellar 
surface are respectively: 100 ms -1 for Model 1; 87ms _1 for 
Model 2; 41 m s" 1 for Model 3; 37 m s" 1 for Model 4. 



where: 



and, 



pr\e r + l{l + m)dr 



2/ + 1 ^-H)! f 1 x - [P r( X )?d X 



2 (f+|m|)W-i 



(19) 



(20) 



£ r and £/j being the radial and horizontal components of the 
eigenfunctions, respectively. In our investigation the functions 
fti(r) have been determined from the numerical integration of 
the Reynolds equation with i max = 4. 

5. Results 

The results of the present investigation are summarized in 
Figs. 0-0 Figure [3] describes, for each of the four com- 
puted Procyon models, the behaviour of the normalized ro- 
tation as a function of fractional radius at different latitudes, 
4> = 0°, 30°, 60°, 90°, for a rotation period P rot = 20 d. The 



normalization of angular velocity is done in terms of a weight- 
ing function based on the moment of inertia: 



n = 



J pr 4 sin 3 efldrde 
f pr 4 sin 3 OdrdO 



with the rotation period P lot = 2n/VL. Figure|4]describes, for 
the same models and rotation rate, the patterns of the merid- 
ional flow in the convection zone. Figures [5] and [6] show the 
same features of the Figs. [3]and|4l but for a rotation period 
-Prot = 33 d. Figure shows the rotational splittings, as de- 
duced from the four differential rotation models of Procyon for 
P rot = 20 d, calculated for the acoustic oscillation modes with 
I = 1, to = 1 and I = 2, to = 1, 2. 

In general, the overall characteristics of distribution of an- 
gular velocity in all the models are of solar type, with an equa- 
torial acceleration and the surface rotation essentially main- 
tained throughout the whole convection zone. 

The total horizontal shear at the surface, (l = (£l oq — 
^poic)/^X ranges from 0.1 to 0.45, depending on models and 
rotation period, that is of the same order of the solar one, 
fi ~ 0.25. The surface shear appears to be stronger in the 
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Fig. 5. The same as in Fig.[3]but for a rotation period P IC 
33d. 



higher rotation case (P ro t = 20 d). For Models 1 and 2, inde- 
pendently of their rotation rate, the rotation latitudinal shear is 
more pronounced at the surface than at the base of the convec- 
tion zone, where there is a tendency to uniform rotation, as in 
the case of the Sun (Figs.[3]and[5j. 

The isorotation surfaces, not presented here, are practically 
the same for all the models and the two rotation rates con- 
sidered, and resemble the solar ones, nearly radial or slightly 
inclined with respect to the equatorial plane, certainly not of 
cylindrical shape, even in the lowest part of the convection 
zone. In most cases the meridional flow is in a single cell ex- 
tending from the poles to the equator, with the surface motion 
directed toward the equator, except for the Models 1 and 2, at 
their higher rotation rate, where an additional small cell ap- 
pears near the equator with the surface motion direct toward 
the pole (Fig. |4}. The maximum speed of meridional motions 
at the surface ranges from 16 ms -1 to 100 ms -1 , depending 
on the models and rotation period (Figs.|4]and|6j. For compari- 
son, the Sun's meridional flow at the surface is directed towards 
the poles with a speed of ~ 10ms -1 at intermediate latitudes, 
as deduced from Doppler shift observations (Hathaway 1996). 

As far as the rotational splittings are concerned, the re- 
sults reflect essentially the different dynamical behaviours of 
the Models 1 and 2, which have deeper convection zones and 
stronger horizontal rotational shears, with respect to those of 
Models 3 and 4, which have shallower convection zones and 
weaker horizontal rotational shears. For the modes with m = 1 






Fig. 6. The same as in Fig.|4]but for a rotation period P rot = 
33 d. Here the meridional motion maximum speeds at the stel- 
lar surface are respectively: 65 m s -1 for Model 1 ; 57 m s -1 for 
Model 2; 17 ms -1 for Model 3; 16 m s -1 for Model 4. 
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Fig. 7. The frequency rotational splittings of acoustic oscilla- 
tion modes, Sv, I = 1, m = 1 (top panel) I — 2, m = 1 
(middle panel) and I = 2, m = 2 (bottom panel), as functions 
of the mode frequency, v, for the four models of differential ro- 
tation with a rotation period P lot — 20 d. The symbols identify 
the models: diamonds: Model 1; squares: Model 2; triangles: 
Model 3; crosses: Model 4. 



(top and middle panels or tig. |7) there is a clear tendency or 
splittings related to the models with weaker rotational shear to 
approach the rigid rotation frequency of 579 nHz, correspond- 
ing to P rot = 20 d, while the splittings related to the mod- 
els with stronger rotational shear are slightly smaller (560 — 
570 nHz), in the case of I = 1, and smaller (520 — 530 nHz), in 
the case of / = 2, than the rigid rotation frequency. Opposite is 
the case of modes with I = 2, m = 2 (bottom panel of Fig.0 
in which the splittings frequencies related to the models with 
weaker rotational shear appear to be smaller than those related 
to models with strong rotational shear. The origin of the dif- 
ferences in the splitting for models 1, 2 and 3, 4 comes mostly 
from the rotational profile, in particular from the latitudinal de- 
pendence of the shear. In fact, if the latitudinal dependence of 
the rotational law is artificially suppressed by setting i max = 
in Eq.(16), the differences in the splittings turn out to be much 
smaller. The fact that a stronger rotational shear lowers the av- 
erage frequency of splittings for modes with m = 1, and vicev- 
ersa enhances the frequency for modes with m — 2 depends 
mainly from the contributions of the angular function terms in 
kernels of Eq.(17), which are negative for m = 1 and positive 
for m = 2 in our models. As a consequence, the sum in Eq.(17) 
is lowered in the first case and enhanced in the second case, this 
trend being more evident for models 1,2 where the latitudinal 
shear is stronger. 

6. Conclusions 

Procyon has already been observed extensively and will pre- 
sumably be observed also in the future with a larger accuracy 
both from ground, by means of multi-site campaigns, or from 
space. Therefore more accurate determinations of the basic pa- 
rameters of the star as well as a higher resolution spectrum of 
its acoustic pulsations are expected, such that can be useful to 
verify the results of the present investigation. 

We have shown that a differential rotation of solar type can 
be maintained even in a very shallow convection zone and that 
its characteristics are enough sensitive to permit the discrimi- 
nation among different stellar models; Models 1 and 2, which 
have a more extended convective envelope, show a larger latitu- 
dinal shear than the Models 3 and 4, which have a less extended 
convective envelope. For the models with the higher rotation 
rate, the meridional motion at the lower latitudes changes from 
poleward, for Models 1 and 2, to equatorward, for Models 3 
and 4. The features of the surfaces of isorotation seem to be ad- 
equate to ensure a possible cyclic dynamo activity in this star, 
whose surface manifestations could depend both on differential 
rotation and meridional flow. 

Precise measurements of the rotational splittings should 
permit in principle the discrimination between models with 
deeper convection zones and strong latitudinal shears (Models 
1 and 2) with respect to those with shallower convection zones 
and weak latitudinal shears (Models 3 and 4). On the basis of 
the present study, a resolution of ~ 10 nHz in the oscillation 
spectrum seems to be necessary to test our models of Procyon. 
Although this resolution is barely reachable for individual os- 
cillation frequency measurements in the near future, even with 
the forthcoming space missions, since the mean value of the 



splitting tor each mode is not strongly dependent on the fre- 
quency, the use of averages over several measured splittings 
could permit to reach the required accuracy. 
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